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Abstract 

KamLAND is a reactor neutrino oscillation experiment with a very long baseline. 
This experiment successfully measured oscillation phenomena of reactor antineutri- 
nos coming mainly from 53 reactors in Japan. In order to extract the results, it is 
necessary to accurately calculate time-dependent antineutrino spectra from all the 
reactors. A simple model of reactor cores and code implementing it were developed 
for this purpose. This paper describes the model of the reactor cores used in the 
KamLAND reactor analysis. 
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1 Introduction 



The KamLAND (A'amioka Liquid Scintillator ^nti-iVeutrino Detector) exper- 
iment [1] is a reactor neutrino oscillation experiment with a very long baseline. 
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The experiment detects low energy antineutrinos (xVs) coming from a number 
of reactors at a typical distance of ~180 km away from the detector. Kam- 
LAND showed evidence for reactor V e disappearance in 2003 [2] , and evidence 
for the distortion of the z/ e energy spectrum, which is consistent with neutrino 
oscillation, in 2005 [3]. In order to extract the oscillation signal, the Kam- 
LAND experiment compares the observed and expected energy spectra. Due 
to the effect of neutrino oscillation, the V e energy spectrum changes as follows 
in the two- neutrino mixing case, 

(Am 2 1 \ 
l-sin 2 20sin 2 ^-J, (1) 

where E Ve is the V e energy, no{E T!e ) is the V e energy spectrum in the absence 
of neutrino oscillation and n(E Va ) is the expected F e energy spectrum with 
neutrino oscillation. Am 2 is the difference of the squared neutrino masses, 
Am 2 = m\ — m\, and L is the distance between the detector and the source. 
The oscillation parameters, Am 2 and sin 2 29, are determined from the fre- 
quency and amplitude of this change. The time dependence of the number of 
z/ e 's per unit power generation changes as much as 10% because the compo- 
nents of nuclear fuel change during the burnup. It is therefore necessary to 
trace the burnup effect for each reactor. However, for the case of KamLAND, 
there are many reactors that contribute V e events and it is practically im- 
possible to calculate the burnup effect using a detailed simulation for all the 
reactors. This paper describes a simple reactor model with which to accurately 
calculate the V e spectrum of each reactor using the routinely recorded reactor 
operation parameters. The parameters include the time-dependent thermal 
output, burnup, and 235 U enrichment of exchanged fuel and its volume ratio. 



2 Antineutrino generation in reactors 

2.1 Expected reactor antineutrino spectrum 

Nuclear reactors are rich V e sources. The fission products are generally neutron 
rich nuclei and undergo successive /3-decays, each yielding one V e . About 6 77 c 's 
are produced per fission along with an average energy release of ~ 200 MeV. A 
typical reactor operating at 3 GW thermal output produces ~ 7 x 10 20 z/ e 's per 
second. At the same time, nuclear transmutation through neutron absorption 
and beta decays of the fuel elements changes the composition of the nuclear 
fuel. In particular, plutonium is created due to neutron capture by 238 U. The 
plutonium production scheme is shown below. The produced 239 Pu and 240 Pu 
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are also fissile isotopes and contribute to the power and u e generation; 
238 U(n, 7 ) 239 U ri/2=23 ' 5min , 239 Np ri/2=2 ' 357d , 239 Pu(n, 7 ) 240 Pu(n, 7 ) 241 Pn. 

v ' E max =1.265MeV B max =0.722MeV v ' v ' 

The four main isotopes, 235 U, 238 U, 239 Pu and 241 Pu, contribute more than 
99.9% of the total power generation. Contributions from other elements can 
be safely ignored. The component of nuclear fuel changes as a function of 
"burnup", which is defined as the time integrated thermal output W(t) per 
initial nuclear fuel mass M, 
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The commonly used unit is giga-wattx day/ton. The burnup is a basic pa- 
rameter which indicates the condition of nuclear fuel. For constant thermal 
output operation, the burnup is proportional to the operation time. Typical 
nuclear fuel resides in the reactor core for 3-5 years and total thermal out- 
put is designed to be 30 MW/t. Therefore, at the end of about one year of 
each operation cycle, incremental burnup from the beginning of cycle (BOC) 
reaches approximately 10 GWd/t. 

The production rate of reactor z/ e 's is obtained from the fission rates and the 
I7 e energy spectra per fission. In the KamLAND data analysis, we compare 
the experimental data with the expected neutrino flux of a superposition of 
all the reactor cores in the integral live time period of the detector. The energy 
spectrum of the expected reactor 17 e flux is written as, 

^(£fJ= E (l-P(Am 2 ,sin 2 2^,E Fo ,L rcactor ))— \ 

reactor " -^reactor 

X X/ VWope(-^F c ) / dt /reactor (t) , (3) 

isotope ^livetime 

where isotope i^h e ) ls the energy spectrum per fission of each fissile isotope, 
-^reactor is the distance from the reactor to KamLAND, P(Am 2 , sin 2 2#, Ep e , -^reactor) 
is the oscillation probability, and /reactor (0 is the time-dependent fission rate 
of each isotope in each reactor. 



2.2 Antineutrino detection 



The KamLAND detector consists of 1 kton of liquid scintillator surrounded 
by 1879 17/20-inch-diameter photomultiplier tubes. In the liquid scintillator, 
I7 e 's are detected through the inverse /3-decay reaction with protons, 

77 c + p -» e + + n. (4) 
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Fig. 1. (a) Reactor I7 e energy spectra for four main fissile isotopes [5]. The shaded 
region for the isotopes gives the uncertainty in the spectrum, (b) Cross section of 
the inverse /3-decay reaction [4]. (c) z7 e observed no-oscillation spectrum for each 
fissile isotope; this is a convolution of (a) and (b). 

The cross section is closely related to the neutron lifetime and known with 0.2% 
accuracy [4]. The energy threshold of this reaction is 1.8 MeV and the z/ e 's 
produced in the /3-decays of 239 U and 239 Np in nuclear fuel do not contribute to 
the reaction. The positron and its annihilation gammas produce scintillation 
light which is proportional to E Vc — 0.8 MeV. The V e visible energy spectrum 
can be written as, 

n(E Vc )=N p a(E Vc )^(E Vc ), (5) 

where N p is the number of target protons in the detector and <j(E Vc ) is the 
inverse /3-decay cross section. Fig. 1 shows the neutrino energy spectrum of 
each fissile isotope, the inverse /3-decay cross section, and the produced visible 
energy spectra in the detector. 





Isotopes 


235yj 238tj 239p u 241p u 


AW< SU 
a 

i 


1 1.52 0.60 0.87 
0.56 0.08 0.30 0.06 
0.61 0.13 0.20 0.06 



Table 1 

Comparisons of the numbers of observed no-oscillation I7 e 's per fission above the 3.4 
MeV V e energy threshold for the isotopes to 235 U, N^ e /N- u . a are the contribu- 
tions of the isotopes to the total number of fission rates in a typical reactor core, 
corresponding to the contributions to the total energy release. £ are the contribu- 
tions of the isotopes to the total number of observed no-oscillation 77 e 's above the 
3.4 MeV energy threshold from a typical reactor core given from Ny c /N v u and a. 
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Above the 3.4 MeV V e energy KamLAND analysis threshold, the comparisons 
of the observed no-oscillation z7 c 's are listed in Table 1. There are significant 
differences in the number of detected 77 e 's per fission. This means that the V e 
spectrum depends on the fuel component and burnup. The number of observed 
events from a typical reactor core decreases by ~ 10% from the BOC to 
b = lOGWd/t. Thus when calculating the V e flux, it is necessary to know 
the dependence of the fuel components on the burnup. The contributions of 
the isotopes to the total number of fission rates in a typical reactor core are 
also listed in Table 1. It must be noted that although 238 U generates a larger 
number of neutrinos per fission, the contribution to the number of observed 
V e events is only ~ 10%. The shaded regions of the neutrino spectra shown 
in Fig. 1 indicate the uncertainty in the spectra. Accordingly, the uncertainty 
in the total number of observed I7 e 's above the 3.4 MeV energy threshold is 
2.5%. 



2. 3 Commercial reactors in Japan 

During the measurement period of KamLAND in Ref. [3] (from 9 March 2002 
to 11 January 2004), 52 commercial reactors in 16 electric power stations and 
a prototype reactor operated in Japan. All Japanese commercial reactors are 
light water reactors (LWRs), 29 are boiling water reactors (BWRs) and 23 are 
pressurized water reactors (PWRs). Both types of LWRs use 3-5% enriched 
uranium fuel. Generally, reactor operation stops once a year for refueling and 
regular maintenance. During the refueling, one fourth of the total nuclear fuel 
is exchanged in BWRs and one third in PWRs. 

To calculate production rates of reactor z/ e 's, knowledge of the correlation be- 
tween the "core thermal output" and the fission rates is required. The core 
thermal output is defined as the thermal energy generated in the reactor cores, 
and it is calculated by measuring the heat balance of the reactor cores. The 
heat taken out by the cooling water, Qf w , is the dominant dissipation source of 
the reactor energy. Other contributions are less than 1%. Therefore, the uncer- 
tainty of the calculated core thermal output is dominated by the accuracy of 
measuring Qf w which itself is dominated by the accuracy of measuring the flow 
of the coolant. The accuracy of the flow of the coolant in turn is determined 
by the uncertainty of the feedwater flowmeters, which are calibrated to within 
2%. In the KamLAND experiment, a value of 2% is used as the uncertainty 
of the core thermal output. 

All Japanese reactors have a contribution of more than ~ 0.1% to the total 
reactor T7 e flux at the location of KamLAND and about half of reactors con- 
tribute between 1% and 7%. Therefore, to accurately calculate the total V e 
flux in the KamLAND experiment, it is required to trace the time variation 
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of the fission rate of all the reactors. To calculate /reactor (t), it is necessary 
to understand the burnup process of nuclear fuel. The process of burnup is 
complicated and depends on the type of core, history of the burnup, initial 
enrichment, fuel exchange history, etc. Detailed simulations exist that calcu- 
late the change of the fuel components in accordance with the burnup. The 
simulation uses the "reactor core analysis method", which traces the burnup 
effect of the three-dimensional fuel component in the reactor core. Ideally, 
it is desirable to perform the simulation for all the reactor cores. However, 
it is practically impossible to perform such a detailed simulation for all the 
commercial reactors, because it is very labor intensive. Based on the above 
circumstances, we have developed a simple core modeling scheme and calcu- 
lation code which can be used to calculate /reactor (t) easily without reducing 
accuracy. We required that the discrepancy of the V e energy spectrum from 
our simplified model be less than 1% from the detailed method. In our model, 
the burnup effect of the fission rate of each fissile isotope in the entire core is 
approximated phenomenologically using reactor operation parameters of the 
nuclear reactor and calculated based on a reference reactor core. All of these 
parameters have been recorded regularly by the nuclear power station, so we 
can use them to calculate the fission rates without requiring additional effort 
from the electric power company. The details of this method are described in 
the next section. 



3 Fission rate calculation 

3.1 Reactor core analysis 

To study the burnup effect of the fission rate in nuclear fuel, sample reactor 
cores under actual operating conditions listed in Table 2 were analyzed with 
the detailed core simulation. These reactor cores were selected to represent 
typical reactor cores in Japan and the target nuclear fuel was chosen to repre- 
sent typical uranium based nuclear fuel including the initial reactor operation 
periods. This analysis used the Core Management System (CMS) codes from 
Studsvik of America, CASMO[6]/SIMULATE[7]. This system performs core 
calculations by combining "two-dimensional multi-group fuel-assembly" anal- 
ysis and "three-dimensional few-group full core" analysis. These codes have 
been extensively compared with measurements. The comparison [8] of the cal- 
culated isotopic concentrations provided by the CMS codes with experiments 
was carried out for the spent nuclear fuel discharged from a BWR type re- 
actor core in Japan. The fuel pin averaged discrepancies between calculated 
and measured isotopic concentrations are less than 7% for the four main fissile 
isotopes. According to the contributions of the fissile isotopes to the number of 
observed z7 e 's from a typical reactor core listed in Table 1, these discrepancies 
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Table 2 

Parameters of new nuclear fuel in six reference reactor cores used in the present 
work. Shown are fuel cycle, average enrichment of new fuel e(%) and volume ratio 
of new fuel V{%). 

of the isotopic concentrations correspond to less than ~ 1% of the number of 
observed z7 e 's per unit energy release. 

The results of calculations of the burnup effect on the relative fission rate of 
each fissile isotope in the BWR type reactors listed in Table 2 are shown in 
Fig. 2. Hereafter, burnup will be defined as the incremental burnup from the 
BOC. From Fig. 2, it is clear that the dependence of burnup can be classified 
in two groups: reactors in their initial fuel cycle, right after commissioning, 
and reactors in their second and higher fuel cycle. Reactor cores in the initial 
operation period, which starts without plutonium, are called "Initial cores". 
In the initial core, the fission rates of 239 Pu increase rapidly from zero after the 
BOC, and fission rates of 241 Pu increase gradually following 239 Pu. After a few 
fuel exchanges, the burnup effect in a particular operation period is indepen- 
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Fig. 2. Relative fission yields for typical reactor cores using a detailed core com- 
position calculation; The ratio is the fission rate divided by the total fission rate 
per unit energy generated. Each line is the burnup dependence of a reactor core, 
multiple cores are shown; the dot-dashed red lines show reactor cores in their initial 
reactor cycle (after commissioning), the dashed blue lines show reactor cores after 
the first refueling (in their second cycle) and the solid black lines show cores from 
the third cycle onwards (equilibrium cores). 

dent of the fuel cycle. This is called "Equilibrium core" . Because the burnup 
effect is large in initial cores and small in equilibrium cores, the dependence of 
the burnup effect on the nuclear fuel parameters is calculated separately for 
the two cases. 



3.2 Initial cores 



The burnup dependence of the number of reactor 77 e 's between 1.0 MeV and 
8.5 MeV V c energy produced in typical initial reactor cores is shown in Fig. 3. 
To extract the dependence on the enrichment, a core with a high enrichment 
of 3.4% is included in this comparison. As shown in Fig. 3, the number of neu- 
trinos decreases with the burnup because of variation of the fuel components. 
In highly enriched fuel, the speed of the decrease is slow. It can be explained 
by the fact that 235 U in highly enriched fuel has a larger contribution to the 
total fission rate. Therefore the enrichment dependence of the burnup effect 
can be represented by transforming the burnup effect in a specific nuclear fuel. 

To study the effect of fuel enrichment, we vary e to first order, substituting 
e o —* e o + Ae, the number of generated z/ e 's can be written as 

N (e , b) = N (e + Ae, b + /3(e , b) Ae) , (6) 
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Fig. 3. (a) Number of reactor z7 e 's in typical initial reactor cores calculated with 
the detailed core simulation. The solid blue line is a core with an enrichment of 
2.2%, the dashed magenta line is 2.5% and the dot-dashed orange line is a high 
enrichment core with an enrichment of 3.4%. (b) Ratio of the number of F e 's in 
two typical initial reactor cores to the number of I7 e 's in a specific core with an 
enrichment of 2.2%. 

where b is the burnup and 

PM= yj?M / 0jm. (7) 

If we know the value of (3(e ,b), the number of z7 e 's can be estimated using 
a function of burnup for a reference value of enrichment, N(eo,b). Here, the 
normalization factor of the burnup, k, is defined as, 

k= 1 + %^Ae. (8) 
b 

Fig. 4 shows values of k in typical reactor cores based on the burnup effect 
in Fig. 3. The factor k is approximately constant with respect to the burnup 
over the complete fuel cycle (Fig. 4(a)) and can be represented by a linear 
function of the new fuel enrichment (Fig. 4(b)). In other words, /3(eo, b)/b can 
be considered to be constant, and the number of TVs is approximated as 

iV(e + Ae,6) = ivL,^. (9) 



In the present work, the value of /3(eo,b)/b is determined by the burnup at 
the end of cycle (&eoc)- The parameters for the actual reactor conditions 
were obtained using the detailed core simulation, k = 1 + 0.65Ae(%) and k = 
l+0.54Ae(%) for the BWR (6 E oc = lOGWd/t) and PWR (6 E oc = 12GWd/t) 
cores, respectively. 
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Fig. 4. (a) Burnup distribution of the k values in two typical initial reactor cores 
plotted as a perturbation of a specific core with an initial enrichment of 2.2%; k is 
approximately constant over the burn cycle, (b) New fuel enrichment dependence 
of the k values. These values are extracted from the results based on the detailed 
core calculation method. 

This simplified reactor core model has so far only been used to calculate the 
number of z7 c 's. However, the neutrino oscillation experiment requires the V e 
energy spectrum. The burnup effect on the V e energy spectrum reflects varia- 
tion of the composition of the fissile isotopes. We assume for the simple core 
model that the fission rates of all fissile isotopes have the same burnup effect 
as when calculating the number of z/ e 's. The model assigns the same correc- 
tions for the new fuel enrichment to the burnup effect to the fission rates. The 
accuracy of this treatment is estimated in a later sub-section by comparing 
the results from our model to the detailed core simulation. The systematic 
uncertainty of the V e energy spectrum is also described there. 



3.3 Equilibrium cores 



The burnup dependence of the number of reactor I7 e 's in typical equilibrium 
cores with exchanged fuel volume ratio of approximately 25% is shown in Fig. 
5. The relative differences of the number of 77 e 's are already less than 1% 
over the complete fuel cycle. To account for the influence of the new fuel, we 
apply the methodology of the fission rate calculation for the initial core to the 
equilibrium core. Nuclear fuel in equilibrium cores contains both new and old 
fuel. The arrangement of new and old nuclear fuel rods is adjusted to burnup 
equally in the reactor core. We treat the fission rate uniformly throughout the 
core. Under the assumption that the fission rates of these two components do 
not have locational dependence, the complete fission rate as a function of the 
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Fig. 5. (a) Number of reactor z7 e 's in typical equilibrium cores calculated with the 
detailed simulation. The solid blue line is a core with an enrichment of 3.0%, the 
dashed magenta line is 3.4% and the dot-dashed orange line is 3.7%. (b) Ratio of the 
number of 17 e 's in two typical equilibrium cores to the number of F e 's in a specific 
core with an enrichment of 3.0%. 

burnup of b in an equilibrium core with the exchanged fuel volume ratio V 
and the new fuel enrichment e is represented as a combination of the fission 
rate of the new fuel and of the rest 



N(b) =N initi&l (b)V + N'(b) 
N'(b) 



V 



l-V 

-^equilibrium (&) — -^initial (&) Vq , 



(10) 



where iV imtia i(6) is the estimated fission rate in the initial core with the fuel 
enrichment of e based on the reference initial core, and -/V equilibrium ((3) is the 
estimated fission rate in the equilibrium core with the exchanged fuel volume 
ratio Vq and the new fuel enrichment e based on the reference equilibrium 
core. Similarly to the initial core calculation, A^ quilibrium (6) is obtained by 
linear approximation of the burnup effect in a reference equilibrium core. The 
values of k are k = 1 + 0.35Ae(%) and « = 1 + 0.29Ae(%) for the BWR 
Oeoc = 10GWd/t) and PWR (6 E oc = 12GWd/t) cores, respectively. 

In addition to the above approximation, accumulation of fissile isotopes pro- 
duced by neutron absorption, particularly 241 Pu produced from the U-Pu 
chain, depends on the averaged absolute burnup of nuclear fuel at the BOC. 
Correction to the fission rate of each isotope is performed as an additional 
contribution using the averaged absolute burnup at the BOC, Absolute 



p/isotope 
J equilibrium I 



( & ) - /equilibrium ( & ) + Af lb0t ° P ° (^absolute) 



where /equilibrium (^) i s the fission rate of each isotope for the incremental burnup 
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from the BOC, b. The correction term is defined as 



A/ isotopc (6 absolutc ) = V isotopc ( b *^ - A , (12) 

\ ^absolute / 

where 71 ^ IS the correction factor for each fissile isotope and ^absolute * s the 
averaged absolute burnup at the BOC in the reference equilibrium core. In 
actual reactor cores, the averaged absolute burnup at the BOC is limited to 
(^absolute /Absolute) ~ 1-4. The effect of this correction to the fission rate is < 5% 
for all fission isotopes except for 241 Pu in BWRs, for which the fission rate is 
corrected by < 30%. 



3.4 Systematic uncertainties 




Fig. 6. Comparison between fission rates in BWR type reactors calculated with 
the detailed core calculation and our simplified calculation. Each line is the burnup 
dependence of a reactor core, multiple cores are shown; the dot-dashed red lines show 
initial cycle reactor cores, the dashed blue lines show second cycle reactor cores and 
the solid black lines show cores from the third cycle onwards (equilibrium cores). 
The largest discrepancy between the detailed core calculation and our simplified 
calculation is for second cycle cores. In practice, this is not of a big concern as 
most of the 52 commercial reactor cores are already equilibrated and the relative 
difference is better than 1%. 

To estimate the accuracy of our simplified calculation method, we compare 
our calculations with the detailed reactor core analysis performed using pa- 
rameters from commercial reactors listed in Table 2. The relative differences 
in the calculated fission rates of each fissile isotopes in BWR type reactors 
are shown in Fig. 6. In this figure, the difference is defined as (/ s Tmpimed(^) ~~ 
/d£3ed)( 6 ) / ^isotope /d2ed( 6 )- Th e difference of the detailed method and sim- 
ple method is less than 3%. At the same time, the relative differences in the V e 
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Fig. 7. Comparison of the v e energy spectrum between the detailed core calcu- 
lation and our simplified calculation in BWR type reactors at the beginning of 
cycle (BOC), middle of cycle (MOC) corresponding to ~ 5GWd/t and end of cycle 
(EOC). Each line is the burnup dependence of a reactor core, multiple cores are 
shown; the dot-dashed red lines show initial cycle reactor cores, the dashed blue 
lines show second cycle reactor cores and the solid black lines show cores from the 
third cycle onwards (equilibrium cores). 



energy spectra for different values of burnup in BWR type reactors are shown 
in Fig. 7. The V e energy spectrum in the equilibrium cores agrees with the 
detailed method within 1%. Because the fission rate of each isotope per unit 
energy generated is balanced between the isotopes, differences in the V e energy 
spectrum are suppressed. The discrepancies are large in the high energy re- 
gion in a initial cycle core and second cycle cores. The maximum discrepancy 
is 2.0% in a second cycle core. However, due to the large number of existing 
Japanese reactors, the V e contribution from second cycle cores is small at the 
KamLAND experiment and can be neglected. The reactor V e flux calculation 
is well performed with an accuracy of 1% in the number of observed z7 e 's in 
all reactor cores, which is smaller than the uncertainty of the V e spectra and 
the reactor core thermal output. 



4 Summary 



To calculate V e flux, reactor cores can be successfully modeled using only a 
few reactor operation parameters. The results of our simplified reactor model 
agree with detailed reactor core simulations within 1% for different reactor 
types and burnup. This error is taken into account in the KamLAND reactor 
neutrino analysis. The simplified model may be applicable to future long- 
baseline reactor neutrino experiments which make use of several reactors. 
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